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ABSTRACT 

We describe a new class of Type I quasars with narrow H/3 broader than 1200 
km/s, above the velocity believed possible for gas in the quasar narrow-line region. 
We identify this group of quasars as a distinct population because of a variety of spec- 
tral and photometric signatures common to these "anomalous narrow- line quasars" 
(ANLs) but atypical of other quasars. One prominent signature is suppression of the 
[OlIl]A4959, 5007 emission-line, in many cases accompanied by broadening similar to 
the H/3 line. We demonstrate that ANLs comprise at least 11% and most likely ap- 
proximately one quarter of the SDSS Type I quasar population at 0.2 < z < 0.8. 
The most striking feature of ANLs is a strong correlation between narrow H/3 width 
and the width of the broad component of the H/3 emission line. Another feature of 
ANLs is a diminished [Oil] line, which might indicate a connection between ANLs and 
the interstellar mediums of their host galaxies, through reduced photoionization or 
star formation. We find that it is difficult to produce ANLs using the current quasar 
standard model. 

Key words: black hole physics — galaxies: evolution — galaxies: nuclei — quasars: 
general — accretion, accretion disks 



1 INTRODUCTION 

Although quasars are luminous enough to be detected at 
large distances, they are too small to be resolved into mul- 
tiple spatial element s except at very short distances (e.g., 
M arcom et al.l(l2006l )) or when lensed by a foreground object 
( Sluse et al.ll2011f ). As a result, our best tool for mapping the 
region around the central supermassive black hole in higher- 
redshift active galactic nuclei has been spectroscopy. Each 
spectral line is associated with a specific ionization potential, 
and the shape of the line corresponds to the peculiar velocity 
profile of the emitting gas, complicated by parameters such 
as densities, abundances, temperatures, and effects of ra- 
diative transfer. There is a long history of studying quasar 
emissi on lines and their correlations (cf. iBoroson fe Green] 
(1992)), much of which has h elped to develop our stan- 
dard model for Type I quasars (|Sulentic et al.ll2000l ; iGaskelll 
2009). In this paper, we describe a new class of quasars with 
distinct spectroscopic features that do not appear to fit eas- 
ily into that standard model. 

The standard picture built fro m existing data c onsists, 
from the black hole outward, of (cf. IPeterson l|2008f )): 



radii away with an exact radius dependent upon the spin 
of the central black hole. The accretion disk and corona are 
believed to emit the quasar continuum. 

• The "broad-line region" (BLR), ~ 0.1 — 1 pc away, from 
which high-velocity gas produces correspondingly broad 
spectral emission lines, with typical FWHM in the 2000- 
20000 km/s range. Prominent emission lines visible in op- 
tical spectra at z > 0.5 include, from ionization poten- 
tial placing them nearest to the central black hole, CIV, 
a broad compo nent of H/3, and Mgl l . In quasar virial 
mass e stimators (IMcLure fc J arvis 2002s iMcLure fc Dunlod 

alJ 120091; 



2004 IVestergaard fc Peterson! [2 006: W ang et 



Onken fc KollmeierlbOO^ ; iRisaliti et al.ll2009l : lRafiee fc Hail 
20101 ). the velocities of BLR gas are assumed to be predomi- 



• The central black hole 

• A hot accretion disk, beginning near the "innermost 
stable circular orbit" (ISCO), typically a few Schwarzchild 



nantly virial in order to use Kepler's Laws to infer the mass 
of the central black hole. However, CIV in particular may 
also be substant ially broadened by r adiation pressure and 
quasar outflows (jMarconi et al.1 [20091 ). 

• The "narrow-line region" (NLR), ~ kpc away, from 
which lower-velocity gas produces correspondingly narrower 
spectral emission lines, with typical FWHM approximately 
500 km/s„ and in some cases higher velocities due to thermal 
or other effects. This FWHM might correspond to typical 
velocities for gas in the interstellar medium being photoion- 
ized by the central black hole. Prominent narrow lines in 
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the optical include the [OIII] fine structure doublet, a nar- 
row component of H/3, and usually [Oil]. 

• A surrounding host galaxy. For quasars, the galactic 
starlight is typically too faint compared to the quasar to 
be detected directly, but galactic spectral lines are often 
present, including both absorption lines and, for many galax- 
ies, a narrow [Oil] emission line believed to be associated 
with star formation. In a typical quasar spectrum, the [Oil] 
emission is likely dominated by photoionization due to a 
non-thermal continuum as opposed to star formation. 

For s ome individual act i ve galactic nuclei, rev erberation 
mapping (|Bentz et al.l 120091 ; IPeterson et al.l 120041 ') has been 
able to confirm the inner portion of this picture, out to the 
H/3 broad emission line. In a time series of spectra for the 
same object, a increase in the continuum luminosity is fol- 
lowed, often hundreds of days later, by a similar flare in CIV 
and then H/3. Assuming the flare propagates outward at the 
speed of light, the delay can be used to infer a radius to BLR 
spectral lines. 

It should be noted, however, that this simple picture 
is merely a broad overview of features common to most 
quasars. Most individual quasars are observed to deviate 
from this model in any of a wide variety of ways, particu- 
larly where high-quality spectra are available and particu- 
larly outside of the broad-line region. 

There have also been hints that for some quasars, 
this simple picture mi ght be more clearly wrong. The 
iBoroson fc Gre en (1992) principal component analysis un- 
covered Eigenvector 1, indicating that for some quasars, 
changes in specific broad spectral lines are correlated with 
changes in specific narrow lines. The population of narrow- 
line Seyfert 1 galaxies contain active nuclei with no tra- 
ditional broad-line component, but with an intermediate- 
velocity component up to ~ 2000 km/s. It is unclear whether 
this intermediate component is a broadened narrow-line 
region, a weak, high-radius broad-line region, or some- 
thing else entirely. Similar ly, some quasar s contain an in- 
termediate H/3 component jllu et al-lkoosl ). Even the [OIII] 
narrow line (A4959, 5007) has been observed to occasion- 
ally be broader tha n 1000 km/s |Marziani et al.1 1 19961: 
iMcIntosh et~ailll999f) . 

The [OIII] fine structure doublet is an especially useful 
indicator for properties of Seyfert galaxi es and other low - 
luminosity active galactic nuclei (AGNs) (lHao et ali l2005). 
[OIII] is a forbidden line that can only be produced in low- 
density gas where there are not enough free electrons to 
allow competing transitions. As such, [OIII] lines are typ- 
ically found around AGN only in the narrow line region, 
as interstellar medium gas is lower-density than gas in the 
broad-line region. Common line-fitting techniques for the 
nearby H/3 line rely upon the assumption that the widths 
of the narrow H/3 component and [OIII] lines are identi- 
cal, in order to disentangle the complicated combination of 
spectral lines surrounding H/3. The [OIII] lines, in turn, are 
typi cally assumed to be narrow based u pon previou s stud - 
ies l|Hao et al.l [20051 ). For example, the IShen et"aH (|2008l ) 
virial mass catalog fit the H/3 broad component along with 
a narrow component linked to the [OIII] line width and at a 
maximum FWHM of 1200 km/s. 

Although fitting each of these narrow lines as one 
Gaussian, with a common width, is currently a stan- 



dard practice, there is ample evidence that the underly- 
ing physics may be more complicated. Forbidden lines at 
higher ionization potentials do show broader velo city pro- 
files then those with lower ionization potentials |Whittlel 
ll985l ; |Pe Robertis fc Osterbrocklll984l ), indicating that the 
narrow- line region may not be monolithic. [OIII], a narrow 
line heavily studied due because it is typically strong and 
available in the optical up t o z ~ 0.8, may be asymmetric 
with multiple components (|Zamanov et al.l |2002| ; iBorosonl 
2005). This might be furthe r evidence of an i ntermediat e- 
line region, or of outflows (cf. lHall et ail |2002| );IHo1 (|20Qgh ). 

In this work, we introduce a set of quasars identified 
via spectral line fitting, characterized by a broadened "nar- 
row" component of H/3 and often also broadened [OIII] (§[2]). 
It turns out these "anomalous narrow-line quasars" (ANLs) 
comprise a distinct population, with a wide variety of prop- 
erties common to all ANLs but atypical of other quasars. 
These properties are discussed in § [31 an d include many of 
the ch aract eristics previousl y reported bv lBoroson fc Greenl 
l|l992h and iHu et all |2008t ) , as well as several properties 
characteristic of narrow-line Seyfert Is, despite our study 
focusing on Type I quasars often containing a 10000 km/s 
broad H/3 component. In addition, we show that ANL spec- 
tra have several additional, previously unreported properties 
characteristic of this new population. We consider whether 
these properties might be an artifact of our fitting routines 
in § U Finally, ANLs seem not to fit easily into the standard 
AGN model, and possible explanations are considered in § [SJ 
as well as issues of nomenclature. 



2 ANOMALOUS NARROW-LINE QUASARS 



We fit all objects in t he IShen et al.l (120111 ) catalog with a 
pr escription similar tolshen et aL I (|200St ). in turn motivated 
bv lMcLure fc Dunlod (|2004l ): 



(i) Using the catalog redshift, we select the spectrum at 
rest wavelengths [4435,4700] and [5100,5535] A and simul- 
taneously fit the sum of a power-law continuum and an 
Fe template (|Bruhweiler fc Vernerl 120081 ) convolved with a 
Gaussian of variable width. 

(ii) With the continuum and iron lines removed, the H/? 
line is fit with two Gaussians and the [OIII] doublet with 
a pair of Gaussians, with the [OIII] doublet is constrained 
to have the 3:1 amplitude ratio physically required by the 
fine structure transitions involved and the two [OIII] widths 
identical to the width of the H/3 narrow component. Typi- 
cal uncertainties are 2-15% in both H/3 velocities, with more 
luminous quasars and broader components typically best de- 
termined. 

(iii) The final result is checked for contamination by 
nearby Hel lines. If Hel contamination is detected, we con- 
struct a n ew line fit with H el removed. This check is not 
part of the Shc n et al.l (|200Sl ) fitting prescription. 

The lShen et alj (|201ll ) catalog includes best fits for objects 
with very poor signal to noise and low-amplitude lines. We 
discard objects with very low H/3 equivalent width and where 
a fit containing no H/3 line had a better X* /DOF than a fit 
containing an H/3 line. Although the Fell width is often sim- 
ilar to the H/3 broad component, this is not required in our 
fitting and for many ANLs, the two have different velocities. 
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Figure 1. Best-fit sum of continuum and iron (green) and total fit (red) to (a) SDSS J000109. 12-004121. 6, a 'standard' quasar, and (b) 
SDSS J000011. 96+000225. 3, an anomalous narrow line quasar. Both H/3 components are also shown. The residual differences between 
spectrum (black) and fit are shown in gray. 
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Because the Fell amplitude is typically far smaller than the 
stronger of narrow H/3 and [OIII] , alternative templates yield 
similar results. A sample end result (the first well-measured 
example sorted by RA) is shown in Fig. Q] 

Contrary to previous work, we find a substantial sample 
of objects with "narrow" H/3 and [OIII] broader than 1200 
km/s using this same technique. In Fig. [TJ right), the quasar 
SDSS J000011.96+000225.3 is shown to exhibit "narrow" 
lines with a FWHM of 1860 km/s. Further, this is not an 
isolated example, but rather part of a large clas s of objects 
with a broadened H/3 narrow component. In the lShen et al.l 
l|201lf) catalog, 30.4% of objects with discernible H/3 broad 
lines contain a narrow H/3 component broader than 1000 
km/s. 

These anomalous narrow-line quasars can be divided 
up into two groups based upon their [OIII] lines. For some 
ANLs, the [OIII] line is suppressed in amplitude, and broad- 
ened comparably to the narrow H/3 component. Several ex- 
amples of these class of ANLs are exhibited in Fig. [2] For 
other ANLs, the [OIII] amplitude is suppressed but the line 
is not equally broadened, as in Fig. [3] For some quasars, it 
is difficult to clearly determine which of these two groups 
is a better description, because ANLs often have broad H/3 
lines, low-amplitude [OIII], and enhanced Fell (§0), leading 
to difficulty fitting the [OIII] line well. We estimate that ap- 
proximately one third of ANLs show broadened H/3 narrow 
and broadened [OIII] of a similar width, and the other two 
thirds show narrower [OIII] than their H/3 narrow compo- 
nent. In both groups, quasars can be clearly identified as 
ANLs because H/3 is a strong, broad line. 

We note that although previous SDSS line catalogs 
|Hao et al.ll2005l ; IShen et al.ll2008l . l201lh have not reported 
a substantial sample of ANLs, several similar objec ts have 
been found in smaller surveys. M arziam et all i| 19961 ) found 
an [OIII]A5007 FWHM over 1100 km/s for 6 of 52 low- 
redshift AGN. In a sample of 32 high-luminosity quasars 
at 2.0 < z < 2.5, higher redshift s than the ANLs re- 
ported in this work, iMcIntosh et al.l |l999) found 12 with 
[OIII]A5007 FWHM =s 1200 km/s, and another six between 
1000 and 1200 km/s, finding a strong correlation b e tween 
[OIII] FWHM and quasar luminosity. iForster et ail (|200ll ) 
point out that the average [OIII] FWHM is 1150 km/s for 
radio-loud and 1160 km/s fo r radio-quiet qua sars in the 
IMcIntosh et all (|l999l ') sample. lAoki et all (|2005l ) present an 
analysis of two quasars with broad [OIII] that is further blue- 
shifted relative to H/3. What we report is that these quasars 
with broadened [OIII] are part of a distinct and larger class 
within the SDSS quasar catalog, and that there is also a sec- 
ond, larger part of that class with broadened H/3 (narrow) 
accompanied by suppressed but unbroadened [OIII]. As de- 
scribed herein, there are a number of other common proper- 
ties (both spectral and photometric) for ANLs, as detailed 
below, that compel us to define a new class of quasars. 



3 ANL PROPERTIES 

In a typical quasar, the broad-line region is at too high of 
a density to allow forbidden transitions such as [OIII], im- 
plying that these quasars with broader [OIII] (ANLs) might 
have atypical structure. Thus, the investigation into ANLs 
begins by determining whether they are atypical in any other 



3000 




Hfi broad component width (km/s) 

Figure 4. Comparison of the narrow-line FWHM and H/3 broad- 
component FWHM for quasars in the SDSS DR7 catalog. Nar- 
rower velocities are accompanied by a wide range of H/3 widths, 
but ANL narrow-line velocities are well-correlated with broad H/3, 
increasing with increasing broad H/3 width. Contours are drawn 
at equally-spaced number densities. 



ways, perhaps accompanied by other distinguishing features 
that might help provide an explanation. 

Indeed, a comparison between the broad H/3 and 
narrow-line FWHM shows two populations (Fig. [4j) , one with 
narrow H/3 and [OIII] apparently uncorrelated with the H/3 
broad-line width and another with broader [OIII] and nar- 
row H/3 well correlated with the broad H/3 component. The 
Type I quasar sample might thus be composed of two dis- 
tinct classes of quasars: (1) quasars with standard narrow 
lines, of the type that has been previously studied and (2) 
ANL quasars in which the physical process responsible for 
anomalous narrow lines also results in broader quasar "broad 
lines". Because low-density gas is required to allow the for- 
bidden transitions producing the narrow [OIII] , the narrow- 
line region is often assumed to lie far from the central black 
hole. A very strong outflow reaching the narrow line region 
could increase the velocity of [OIII] and the narrow com- 
ponent of H/3, and would at the same time also increase 
the velocity of the broad H/3 component coming from the 
broad-line region closer to the central black hole. 

These two classes may be more difficult to distinguish 
than previously thought. As indicated in Fig. 2] a cut at 
~ 1200 km/s is sufficient to remove standard quasars from 
an ANL population, but objects with low FWHM of both 
H/3 components, classed for our study as standard quasars, 
might instead lie on a continuation of the ANL branch, or 
might be a combination of some objects with standard and 
some objects with ANL physics in their broad-line regions. 
If so, the true ANL population would comprise more than 
30% of the S DSS catalog, and s ome ANLs would have been 
well- fit in the lShen et al.l (|201ll ) catalog. This possibility is 
discussed further in § 15.21 

For further investigation, we divide the SDSS DR7 
sample into bins by narrow-line width, as in Table [3] 
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Figure 2. Four examples of anomalous narrow-line quasars exhibiting [OIII] broadening comparable to the narrow H/3 component. The 
total fit (red) is composed of a continuum plus Fell fit (blue), two H/3 components (green), and [OIII]. (a) SDSS Jlf 2048.01+631453.9, 
z=0.514; (b) SDSS J145142.83+610747.8, z=0.296; (c) SDSS J081528. 88+333254.2, z=0.482; (d) SDSS J095912.92+445059.0, z=0.407 



Table 1. Properties of quasars in the SDSS DR7 catalog binned by narrow-line width, including the fraction showing broad abso rption 
in Mgll (Shen et al. 2011) and the best-fit Fell amplitude (arbitrary units). 



a (km/s) 


FWHM 


Color 


N 


Frac. 


'Best' N 


Frac. 


BAL 


Fell (arb.) 


log Lbol 


100-200 


235-471 


Black 


3549 


0.217 


1162 


0.217 


0.0033 


1.59 


45.41 


200-300 


471-706 


Black 


4941 


0.303 


1018 


0.274 


0.0050 


2.14 


45.51 


300-400 


706-942 


Red 


2309 


0.141 


657 


0.123 


0.0098 


3.25 


45.61 


400-500 


942-1177 


Red 


1366 


0.084 


420 


0.078 


0.0069 


4.14 


45.68 


500-600 


1177-1413 


Yellow 


1086 


0.067 


394 


0.074 


0.0084 


7.10 


45.72 


600-700 


1413-1648 


Green 


1018 


0.062 


432 


0.081 


0.0065 


9.37 


45.82 


700-800 


1648-1884 


Cyan 


743 


0.045 


321 


0.060 


0.0116 


10.48 


45.90 


800-900 


1884-2119 


Blue 


424 


0.026 


207 


0.039 


0.0249 


11.96 


46.00 


900-1000 


2119-2355 


Magenta 


227 


0.014 


110 


0.021 


0.0284 


13.75 


46.01 



Quasars with standard narrow lines are most common, with 
the population peaking near the expected 500 km/s line 
FWHM. However, 30.4% of quasars have narrow-line ve- 
locities greater than 1000 km/s, and in nearly one quarte r 
of quasars these velocities are above the Sh en et all (|2008l ) 
limit of 1200 km/s. Although Fig. [4] indicates that there are 
two populations of quasars, the quasar fraction decreases 
monotonically at velocities above the peak. This may indi- 
cate that quasars are not necessarily either ANLs or "stan- 
dard" quasars, but rather can also exist in an intermediate 
state. Anomalous narrow lines are accompanied, on average, 



by an increased luminosity. However, ANLs are also increas- 
ingly likely to co ntain broad absor ption lines as measured 
from Mgll in the lShen et~al] l|201ll ) catalog. 

For a variety of reasons, our best fit may not be a good 
match for even a well-measured SDSS spectrum, and this 
is one of the limitations of any automated line-fitting pre- 
scription with a small number of free parameters. For well- 
measured spectra where narrow H/3 and [OIII] have similar, 
Gaussian profiles, our fit is generally good. However, many 
quasars best-fit with narrower "narrow" lines have non- 
Gaussian profiles, and many anomalous narrow-line quasars 
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Figure 3. Two examples of anomalous narrow-line quasars exhibiting [OIII] suppression, but without [OIII] broadening comparable to 
the narrow H/3 component. The total fit (red) is composed of a continuum plus Fell fit (blue), two H/3 components (green), and [OIII]. 
(a) SDSS J170208.34+330700.6, z=0.658; (b) SDSS J110103.76+411317.1, z=0.455 
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Figure 5. Comparison of co-added quasar spectra around the H/3 
and [OIII] lines binned and colored as in Table [3] The continuum 
and Fell lines have been subtracted. 



show a mismatch between the H/3 and [OIII] profiles. In such 
cases, our fitting routine will typically be dominated by the 
higher- amplitude of narrow H/3 and [OIII], which typically 
is [OIII] for lower velocities and H/3 at higher narrow-line ve- 
locities (Fig.0 further described below). Some mismatched 
ANLs may be objects with a broad blue [OIII] wing com- 
ponent tied to H/3 and in addition to t he OIII narrow com - 
ponent, similar to those discussed in I Zhang et ail l|201lh . 
Quasars with ~ 2000 km/s narrow H/3 and narrower [OIII] 
are certainly anomalous, but might have a different cause 
than ANLs with ~ 2000 km/s FWHM for both H/3 and 
[OIII]. 

We produce a 'best' sample of quasars for which (1) 
the overall fit x 2 /DOF is good and (2) the x 2 /DOF for the 
~ 50 — 150 pixels within 2a of the narrow H/3 and [OIII] is 
within 0.5 of the total x 2 /DOF. As in Table this sample 
has a slightly larger ANL fraction than the overall catalog. 
A very conservative lower bound on the true ANL fraction 
would be 11%, the fraction of quasars in the entire sample 
that are ANLs in the 'best' sample. A much more likely 
explanation is that ~ 11% of quasars are ANLs with both 
H/3 and [OIII] broadened (Fig. [5]), while an additional ~ 
20% of quasars are ANLs with H/3 broadened but [OIII] 
suppressed and unbroadened (Fig. [3j . 

Spectra within each bin from Table are co-added to 
examine the average dependence of spectral lines on OIII 
width. Each individual spectrum is first smoothed to the 
resolution of the SDSS spectrograph, then normalized to the 
monochromatic 5100 A flux and averaged. Subtracting the 
best-fit continuum and Fell template results in the spectral 
lines shown in Fig. [S] As indicated by Fig. 3] varying widths 
for standard "narrow" lines (black, red) in Fig. [5] are not 
associated with strong changes in the broad component of 
H/3. However, broader "narrow" lines are associated with an 
increase in width of H/3 as well as increasing H/3 line flux. 
When the [OIII] width increases, it is associated with a de- 



clining amplitude, and in total a nearly-constant equivalent 
width. 

Several other spectral features are also responsive to 
increased narrow-line width (Fig. |SJ . The continuum tilt in- 
creases for ANLs, making them bluer than other quasars. 
This bluer tilt occurs despite a higher BAL fraction among 
ANLs (Table which indicates ANLs might be on av- 
erage dustier than other quasars. Perhaps this continuum 
slope is explained by ANLs having a higher tempera- 
ture, also consistent their higher average luminosity. ANLs 
also s how stronger Fell lines (see also iBoroson fc Greer] 
(1992), Table |3j). An effort to find ANL indicators avail- 
able at other redshifts where [OIII] is unavailable finds 
that [NeV] broadens along wi th [OIII], while SII does not 
l|Steinhardt fc Silverman! 120111 ). The stronger and broader 
H/3 line might be evidence that the nature of the quasar 
broad-line region is changing, whether due to a strong out- 
flow or other accretion physics. However, the Mgll line (Fig. 
does not increase in width in ANLs. Since Mgll has an 
ionization potenti al placing it at a larger radius from the 
central black hole |Peterso3[2008), one possible explanation 
might be an outflow propagating only partway through the 
broad-line region. 

However, the [Oil] A3727 A line appears sensitive to 
these changes (Fig. [6}. [Oil] originates in both narrow- line 
region g as around the quasar and s tar formation in the host 
galaxy i|Hol 120051 ; Ikim et alj |2003 ). If broadened [OIII] is 
evidence of an outflow broadening the entire narrow-line re- 
gion, then any [Oil] coming from near the quasar would also 
be broadened. The [Oil] amplitude is declining in Fig. H3 but 
the line width is not increasing. 

The decline in [Oil] amplitude indicates that the com- 
bination of [Oil] from all sources, including star formation, 
diminishes with increasing [OIII] width. One intriguing possi- 
bility is that star formation in ANL hosts is quenched by outflows 
from the central black hole. It is also possible that all quasars in- 
hibit star formation, and that the reduction in [Oil] flux between 
standard Type I and ANL quasars is due to a reduction in the 
photoionized [Oil] in the narrow-line region. 

The [Oil] ccntroid appears to shift between different bins in 
Fig. [6] Although this could be interpreted as a skew in [Oil], the 
redshift is determined primarily by the highest-equivalent width 
lines. If H/3 is skewed (Fig. |5j in ANLs, the result may be a sys- 
tematically incorrect redshift, resulting in locating low-amplitude 
lines such as [Oil] at the wrong rest wavelength. A possible in- 
dication of this skew is the offset between the broad and narrow 
components in H/3 fits to ANLs /Figs. [2l3t . If the narrow compo- 
nent is blueshifted while the broad component lies at the velocity 
of the host galaxy, it might induce this sort of skew. This ef- 
fect may make high-precision redshift determination difficult for 
ANLs, as their broad emission lines are skewed and their galactic 
lines are difficult to measure due to their small amplitude. 

If this redshift error leading to an [Oil] offset is systematic 
and uniform for fixed narrow-line width, the average decrease in 
[Oil] amplitude will be well-indicated by the co-added spectrum. 
However, if the offset is non-uniform, high-amplitude lines with 
different centroids may be co-added to produce a low-amplitude, 
broader line. Therefore, we fit the [Oil] line directly and inves- 
tigate its properties. As with our H/3 and [OIII] fits, we first fit 
a combination of continuum and broadened Fell template, then 
find the best-fit Gaussian for the leftover line profile. We find 
that the [Oil] equivalent width declines by a factor of 2.25 for 
ANLs compared to standard quasars. The best-fit [Oil] width 
remains typical of a narrow galactic emission line, and the de- 
creased line flux comes from to a decrease in [Oil] amplitude. For 
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Figure 6. Comparison of co-added quasar spectra at different narrow-line widths binned and colored as in Table [3] For shorter wave- 
lengths, only spectra at sufficiently high redshift were co-added. At top left, the Mgll region. At top right, the [Oil] region. 



low-amplitude [Oil], noise is often incorporated into the best-fit 
line profile, so the true decline may be larger. 



ing a techn i que p reviously used bv lMcLure fc Dunlod l|2004h and 
Shcn ct al. (2008), we hope to have avoided introducing any novel 
line-fitting errors. However, we should highlight two specific dan- 
gers: 



4 DEPENDENCE OF ANLS ON FITTING 
TECHNIQUES 

Because of the long history of studying quasar spectra, a report 
of a large, new spectroscopic class of quasars discovered via spec- 
tral line fitting should be met with appropriate skepticism, and 
it is important to consider whether ANLs might just be an arti- 
fact of fitting techniques. Ultimately, the H/3-[OIII] complex arises 
from complicated physics but is fit with a small number of pa- 
rameters in an automated way, without the underlying model 
being adjusted individually for tens of thousands of quasars in 
the SDSS DR7 and future catalogs. As a result, most fits should 
contain errors. In an effort to more fully disclose some of these 
errors, we have made a large sample of our fits for ANLs available 
at http://rnernber.ipmu.jp/charles.steinhardt/ANL/ By choos- 



4.1 Continuum Fitting 

Errors in continuum fitting could result in incorrectly determin- 
ing the properties of broad spectral lines, either missing the broad 
component of H/3 if too much flux is considered part of the con- 
tinuum or introducing a false broad component by under-fitting 
the continuum. Although for ~ 500 — 1000 A windows the con- 
tinuum is well-fit by a power law, the slope of that power law 
changes betwee n the Mgll a nd H/3 region in composite spec- 
tra jVanden Berk et al.|[200lh . so it must be locally determined. 
However, most of the pixels near H/3 are contaminated by ei- 
ther Fell or strong spectral emission lines. Our fitting tech- 
nique solves this problem by fitting simultaneously for Fell emis- 
sion and the continuum. An alternative technique would involve 
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choosin g a window free o f Fell contamination, such as 5080- 
5100 A lForster et alj|200ll) , and using these windows to fit the 
continuum alone. However, in an ANL with broadened [OIII], 
5080 Amay not be free of [OIII] contamination, so this window 
could also provide an erroneous continuum slope. 

By simultaneously fitting Fell and continuum, we fit hun- 
dreds of SDSS pixels with four parameters, resulting in a well- 
determined fit. However, although the continuum fits exhibited 
in Figs . ITI31 appear reasonable, it is difficult to prove that any in- 
dividual spectrum has been well-fit. Fortunately, it is the broad, 
rather than the narrow, component of H/3 that will be strongly 
affected by errors in continuum fitting. Our report of a large pop- 
ulation of ANLs relics on fitting of the narrow component, and 
thus should not be susceptible to errors in continuum fitting. 

4.2 Iron Contamination 

Although Fell is a larger problem when fitting Mgll rather than 
H/3, there are also Fell lines extending near [OIII]. Further, the 
strength of iron lines increases strongly in ANLs (Tabic |3). As 
with continuum fitting, Fell contamination should not be a prob- 
lem when determining whether an object is some type of ANL, 
as the narrow H/3 component is stronger than Fell and well- 
measured. However, the combination of strongly suppressed [OIII] 
and augmented Fell emission is more likely to result in errors fit- 
ting [OIILU5007. 

The presence of stronger iron lines does make Fell eas- 
ier to detect and fit. Further, although there are several dif- 
ferent Fell templates available, they are typically derived from 
the spectra of objects with narrower Fell lines, such as 1 Zw 1 
llBruhweiler fe Vc rncr 2008). Because these templates are then 
convolved with a broad Gaussian to simulate Fell with a FWHM 
of thousands of km/s, differences between templates are mini- 
mized, so alternative templates produce a similar ANL sample. 

The reduction from the full ANL sample to the 'best' sam- 
ple presented in § [3] uncovered a large population of ANLs in 
which [OIII] was suppressed but not broadened. Perhaps what 
we have fit as suppressed and broadened [OIII] in these other 
cases is instead dominated by strong Fell emission at broad-line 
region velocities. Because several ANLs show not just a broadened 
[OIII]A5007 line but also equally-broad [OIILU4959, for some ob- 
jects [OIII] is unambiguously broadened (Fig. [2^,,c). In other ob- 
jects, a combination of Hel contamination and highly-suppressed 
[OIIL1A4959 might make an alternative explanation involving Fell 
more plausible (Fig. f2]d). 

The Fell FWHM is often linked to the broad H/3 compo- 
nent FWHM. However, in ANLs with broadened [OIII], the [OIII] 
FWHM is instead linked to the H/3 narrow component, an appar- 
ent mismatch. However, ANL broad-line regions appear to be- 
have anomalously in several ways compared to BLRs in standard 
quasars, so perhaps iron emission is also different. Iron contam- 
ination is therefore a plausible concern when fitting [OIII] for 
some ANLs, but cannot be responsible for miscatcgoring objects 
as ANLs based upon the H/3 narrow component. 

4.3 Modeling H/3 With Two Gaussians 

In this work, we choose to model the H/3 line as a pair of Gaus- 
sians, and each of the two [OIII] lines as one Gaussia n. This is in 
keeping with s tandard practice fo r large c atalogs fcf. IShen et al.l 
ll201ll , 1200^) ; iMcLure fc Dunlod i2004\ : iRafiee fc Halll J2010T) . 
etc.). Even some past work using more complicated fi tting pre- 
scriptions reported their results using two Gaussians l lHu et al.l 
2008). As such, if the correct interpretation of this result turns 
out to be merely that the two-Gaussian fit is faulty for ~ 30% 
of quasars, it would still be important to understand and correct 
for this problem, and that correction would change the results of 



existing statistical studies of the SDSS quasar population using 
H/3 line fitting, including virial mass catalogs. 

However, the question remains as to whether two Gaussian 
components is a good representation of the H/3 line, particularly 
for ANLs. When only one component of an ANL is allowed to 
be broader than 1200 km/s, effectively the entire line shape is 
modeled as one Gaussian, where the fitting in this work uses 
two. Fitting two broad Gaussians to a line shape closer to one 
Gaussian, but not exactly Gaussian, might lead to fit composed 
of the dominant Gaussian and a second, correlated Gaussian that 
would be artificial and representative of line non-Gaussianity. So, 
is one Gaussian actually a better fit for ANLs? For th at matter, 
is a th ird Gaussian an even better fit, as reported by iHu et all 
ll2008h ? 

The simplest measure of whether two Gaussians provide a 
better fit is to examine the x 2 /DOF for each fit. We consider the 
part of the spectrum within 2cr of the peak, as defined by the 
broadest component. For all quasars, the x 2 /DOF averages 1.31 
for two-Gaussian fits and 1.70 for one-Gaussian fits. For ANLs, 
two-Gaussian fits average 1.30 while one-Gaussian fits average 
2.02. Adding a third Gaussian, however, provides no such sub- 
stantial help; a three Gaussian fit reduces the x 2 /DOF from 1.31 
to 1.29 for all quasars and, like the two-Gaussian fit, averages 
1.30 for ANLs. Thus, it is clear that the second-strongest compo- 
nent of the H/3 lincshape is strong, while the third-largest is much 
weaker. A two-Gaussian fit therefore seems reasonable, although 
the correct line shape may not contain Gaussians at all; a cor- 
rect fit should produce a x 2 /DOF of 1. Howeve r, errors in SDSS 
are li kely correlated between neighboring pixels ( McDonal d et all 
2006), so the reported uncertainties used in producing these value 
for x 2 may not be trustworthy. This x 2 analysis is a good sanity 
check on our fitting technique, but is certainly not proof that the 
lines are being fit well. 

Ultimately, fitting lines is an art, there are many com- 
plex prescriptions that can be used to produce higher-quality 
fits for individual spectra, and it is unrealistic to tailor indi- 
vidual fits for 100,000 spectra in SDSS, or even more in fu- 
ture surveys. We believe the best we can do as scientists is 
to use public data, make a collection our fits public rather 
than only exhibit our best examples (a collection is available at 
http://member.ipmu.jp/charles.steinhardt/ANL/), and encour- 
age other experts to use alternative prescriptions for H/3. It ap- 
pears very difficult to fit ANLs well with one broad component 
and one component under 1200 km/s in FWHM. 

In general, there are many errors arising from fitting the 
spectrum of a complex object with a simplistic model containing 
little physics and a small number of parameters. These errors will 
likely result in some incorrect fits, it would be difficult to believe 
reported H/3 narrow widths to within 1%. However, many ANLs 
have narrow H/3 measured at over 2000 km/s, and even a rudi- 
mentary fit "by eye" of the narrow H/3 lines displayed in Figs. 
IH3I cannot accommodate typical narrow H/3 velocities of under 
1000 km/s. These details of line-fitting are important when, for 
example, using precision line measurements to produce a quasar 
virial mass. However, they cannot be responsible for overestimat- 
ing the narrow H/3 width by such a large factor in 30% of the 
SDSS quasar catalog. 



5 DISCUSSION 

We report a new class of Type I quasars with broadened nar- 
row H/3 and [OIII]. These quasars (ANLs) include over 1/4 of all 
~ 16000 SDSS quasars for which H/3 is well-measured (z < 0.8). 
However, previous line-fitting catalogs of SDSS quasars have not 
reported finding a large population of ANLs, likely because se- 
lection and fitting routines assume that [OIII] is a narrow line. 
The existence of these objects now necessitates a re-examination 
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of other selection criteria in SDSS and similar surveys. In partic- 
ular, the presence of a strong, narrow [OIII] line is often a strong 
component of Type II quasar and Seyfert galaxy selection, since 
broader lines may not be present. Any type II ANLs or ANL 
Seyferts would likely be selected against during the production of 
such catalogs. 

Several previous examples have been reported of other AGNs 
with anomalous narrow line s. These include S eyfert s with broad- 
ened or skewed [OIII] (cf. IXu fc Komossal <201ll)1. as well as 
quasars with multiple broad H/3 components jHu et alJbOQSft . 
with the broa der broad componen t redshifted with r egard to the 
narrower one JSulentic et al.ll2002l; ICollin et ai1l200rj) . 

ANLs, however, not only show two well-centered broad com- 
ponents but represent approximately one quarter of the SDSS 
quasar catalog at z < 0.8, with thousands of examples. As such 
their origin is a compelling puzzle. The spectra of these quasars 
provide enough evidence to analyze several possible explanations: 
Could ANLs be created by an outflow into the narrow- 
line reg ion? Given that outflows are a common property of most 
quasars ( Elvis 2000J) , it would not be surprising to find that a large 
population exists with discernable evidence for possible non-virial 
gas motion in the broad-line (and narrow-line) region. Although 
there is a slight skew to [OIII] and the narrow component of H/3 
in co-added ANL spectra (Fig.O, the lines are mostly symmetric. 
Thus, any outflow must be similarly symmetrical, meaning that 
light from both sides of the outflow is visible. Otherwise, these 
lines would be more strongly ske wed, just as stro ng asymmetric 
outflows can skew CIV emission (Shen et alj|2008h . 

Moreover, Mgll, which emanates from the outer portion of 
the broad-line region, does not increase in width in ANLs (Fig. 
l6|l. The narrow- line region, containing [OIII], lies further from 
the central black hole than the broad-line region. Thus, since any 
outflow does not reach past the edge of the broad-line region, 
ANLs are not created by an outflow into the narrow-line region. 
The broad component of H/3 also comes from the broad-line re- 
gion, but its ionization potential places it closer to the central 
black hole than Mgll. Therefore, an outflow that does not reach 
Mgll-emitting gas could still be responsible for higher velocities in 
gas producing the broad component of H/3. An alternative option 
might be an outflow, larger than the accretion disk, that is able 
to emit H/3 (and, presumably, CIV), but unable to emit Mgll. 
Could ANLs be created by [OIII] emission in the broad- 
line region? As a forbidden line, [OIII] requires a low density. 
In a typical Type I quasar, the broad-line region is too dense 
to allow [OIII] emission, which is why [OIII] does not have a 
prominent broad-line component. In ANLs, the line flux from H/3 
and Mgll, both emitted in the broad-line region is comparable 
to or larger than other Type I quasars. So, the broad-line region 
density is unlikely to be substantially lower. Further, if [OIII] and 
the narrow component of H/3 are in the broad-line region, where 
is the broad component of H/3 located? 

One way to test this picture is to co nsider the evidence from 
virial mass estimates. As described in l|Steinhardt &: Silverman! 
2011), the H/3 broad component yields a systematically larger 
mass estimate than Mgll for ANLs. Unlike H/3, Mgll is not 
strongly correlated with changes in [OIII] for ANLs, and there- 
fore Mgll is likely the better mass indicator. For agreement, the 
H/3 broad component must be placed at less than half of the ra- 
dius inferred from the continuum luminosity - broad line region 
radius relation calibrated by reverberation mapping. If the H/3 
narrow component and [OIII] were instead placed at that radius, 
it would result in virial masses 1-2 dex below those produced us- 
ing Mgll. Thus, [OIII] emission from virialized broad-line region 
gas does not seem to fit the available evidence, although some of 
that evidence may not be reliable. 

Could ANLs be typical Type I quasars viewed only from 
certain angles where an outflow is seen most clearly? 

[Oil] emission comes in part from star formation and in part from 



larger radii than [OIII] , the latter due to its lower ionization state. 
Thus, [Oil] emission sharply declining (Fig. |6{ with broadened 
[OIII] seems to indicate that ANLs lie in host galaxies, or at least 
~ kpc scale central regions, with different properties than other 
Type I quasars. This would not happen merely due to geometry. 
Indeed, it is a puzzle that the [Oil] line indicates a connection 
between the host galaxy and the quasar, while the unchanging 
Mgll indicates that the connection is not directly due to a direct 
outflow. 

Could ANLs be a pre-turnoff phase of the quasar duty 
cycle? The decreased [Oil] equivalent width, coming from a com- 
bination of star formation and narrow-line region gas, might in- 
dicate a reduced supply of gas and dust to the central black hole, 
perhaps further reduced by a strong outflow. This, in turn, could 
cause the central quasar to go quiescent, with ANLs an interme- 
diate phase as the quasar turns off. However, the average lumi- 
nosity of ANLs is actually higher than other quasars (Table [3| ■ 
An upcoming luminosity decrease would not necessarily render 
the quasar quiescent, and the response time to a reduced inflow 
~ kpc away is long. Nearby, ANLs might actually have more dust 
available than average Type I quasars, as they are more likely to 
contain broad absorption lines (Table l3t . This picture also does 
not yet explain the broadened H/3 and [OIII] lines. 

So, although several scenarios present themselves, none yet 
seems able to fit all of the available evidence. 



5.1 ANLs and Feedback 

There is known to be a link between the properties 
of lower-luminosity active galactic nuclei and their hosts 
teauffmann et alj|2003t ISilverman et aljlioblil : ISchawinski et al.l 
l20ld ; ICardamone et all |2010T) . Although co-evolution of the 
quasar and its host galaxy has long been expected and even as- 
sumed, direct evidence for this idea has not yet been forthcoming, 
primarily because the quasar dominates its host and because the 
peak quasar number density lies at z ~ 2, where resolving the 
host galaxy as separate from the quasar cannot be done with cur- 
rent ground-based instruments and where the quasar itself, the 
accretion disk and the broad-line region cannot be resolved with 
any existing instruments. 

One particularly striking feature of ANL quasars is the 
strong link in co-added ANL spectra (Fig. [6j between [OIII], pro- 
duced either in nearby heated interstellar medium or closer to 
the central black hole, and [Oil], produced in photoionization ~ 
kpc from the black hole and in star formation in the host galaxy. 
The correlation between ANLs and their hosts appears to be evi- 
dence of some sort of feedback or linked evolution. Further work, 
including a physical model for that link, is needed to understand 
whether we are seeing feedback, co-evolution, perhaps some other 
process that, as a byproduct, both quenches star formation and 
broadens the [OIII] line. 

5.2 Nomenclature 

One of the difficulties in reporting anomalous narrow-line quasars 
has come in choosing the terminology. In many ways, a bet- 
ter description might invoke an unsatisfying combination such 
as "broad narrow-line", or even "broad narrow H/3", in describ- 
ing these quasars. Alternatively, perhaps they might be named 
in analogy with Seyfert nomenclature (e.g., NLSyl). It was even 
pointed out to the authors that with ANLs comprising as much 
as 30% of the quasar catalog, their narrow lines aren't really 
anomalous. The root problem is that a set of lines convention- 
ally referred to as "narrow" are often found at velocities typically 
termed "broad" , and it may be difficult to satisfactorily describe 
these objects without altering common terminology. It is our hope 
that future work will produce a physical understanding of ANLs 
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and allow a name representative of their origin rather than their 
spectra. 

It should also be noted that the selection of ANLs in this 
study has been based entirely upon spectroscopic criteria, i.e., 
their H/9 narrow component is anomalously broad. A more nat- 
ural categorization might be based upon the physics of their 
broad-line regions. For example, quasars might be classed as ei- 
ther virial (VQ) or non-virial (NVQ) based upon their broad-line 
regions. As discussed in followu p work, many ANLs appear t o 
have non-virial broad-line regions jSteinhardt fc Silvermanl201ll) . 
One might guess that quasars on the diagonal branch of Fig. [4] 
might be NVQs, while quasars on the horizontal branch are VQs, 
with the quasars in the lower-left intersection of the two branches 
perhaps consisting of a mixture of both VQs and NVQs. In this 
interpretation, what we have currently labeled our ANL sample 
would consist most of NVQs with a few VQs mixed in, while our 
"standard" quasar sample would consist mostly of VQs but with 
a subsantial number of NVQs as well at lower-left. 

This work has defined ANLs as a population for which stan- 
dard assumptions about the broad- and narrow-line regions ap- 
pear to be false. This is a useful definition because line-fitting 
catalogs are often predicated upon these assumptions, and for 
ANLs an amended technique must be used. The best definition 
would simultaneously divide quasars spectroscopically and physi- 
cally. However, if additional spectral properties cannot break the 
degeneracy in the lower-left of Fig. [4] it may even be possible that 
multiple categorizations are required. 
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